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Big-bang nucleosynthesis is one of the cornerstones of the standard cosmology. 
For almost thirty years its predictions have been used to test the big-bang model 
to within a fraction of a second of the bang. The concordance that exists between 
the predicted and observed abundances of D, 3 He, 4 He and 7 Li provides impor- 
P., 1 tant confirmation of the standard cosmology and leads to the most accurate de- 

termination of the baryon density, between 1.7 x 10 31 gcm 3 and 4.1 x 10~ 31 gcm -3 
(corresponding to between about 1% and 14% of critical density). This mea- 
surement of the density of ordinary matter is crucial to almost every aspect of 
cosmology and is pivotal to the establishment of two dark-matter problems: (i) 
most of the baryons are dark, and (ii) if the total mass density is greater than 
about 14% of the critical density as many determinations now indicate, the bulk 
of the dark matter must be "nonbaryonic," comprised of elementary particles left 
from the earliest moments. We critically review the present status of primordial 
nucleosynthesis and discuss future prospects. 



1 Introduction 

Because of the extremely high temperatures that existed during the earliest moments it was 
too hot for nuclei to exist. At around 1 sec the temperature of the Universe cooled to 10 10 K, 
and a sequence of events began that led to the synthesis of the light elements D, 3 He, 4 He 
and 7 Li. The successful predictions of big-bang nucleosynthesis provide the earliest and most 
stringent test of the big-bang model, and together with the expansion of the Universe and the 
2.726 K black-body cosmic background radiation (CBR) are the fundamental observational 
basis for the standard cosmology. 

Big-bang nucleosynthesis began with the pioneering work of Gamow, Alpher, and Herman 
who believed that all the elements in the periodic table could be produced in the big bang 
|75| ; however, it was soon realized that the lack of stable nuclei of mass 5 and 8 and Coulomb 
repulsion between highly charged nuclei prevent significant nucleosynthesis beyond 7 Li. In 



1964, shortly before the discovery of the CBR, Hoyle and Tayler |76| argued that the big 
bang must produce a large 4 He abundance (about 25% by mass) and this could provide the 
explanation for the high 4 He abundance observed in many primitive objects. At about the 
same time, Zel'dovich realized that a lower temperature for the Universe today implied a 
greater mass fraction of 4 He produced in the big-bang, and concluded that the big-bang model 
was in trouble. While his reasoning was correct, through a comedy of misunderstandings he 
mistakenly believed that the primeval mass fraction of 4 He was at most 10% and that the 



temperature of the Universe was less than about 1 K |77|. 

After the discovery of the CBR by Penzias and Wilson in 1965, detailed calculations 
of big-bang nucleosynthesis were carried out by Peebles |7J| and by Wagoner, Fowler and 



Hoyle |fr9| . It soon became clear that, as Hoyle and Tayler had speculated, the origin of the 
large primeval fraction of 4 He was indeed the big-bang, and further, that other light elements 
were also produced. However, the prevailing wisdom was that D and 7 Li were produced 
primarily during the T Tauri phase of stellar evolution and so were of no cosmological 
significance |80| . The amount of 4 He produced in the big bang is very insensitive to the 



cosmic baryon — that is, ordinary matter — density (see Fig. 1), and thus it was not possible 
to reach any conclusions regarding the mean density of ordinary matter. 

Since the other light elements are produced in much smaller quantities, ranging from 
10" 5 or so for D and 3 He to 10~ 10 for 7 Li (see Fig. 1), establishing their big-bang origin 
was a much more difficult task. Further, it is complicated by the fact that the material we 
see today has been subjected to more than 10 Gyr of astrophysical processing, the details 
of which are still not understood completely. However, over the past 25 years the big-bang 
origin of D, 3 He, and 7 Li has been established, not only further testing the model, but also 
enabling an accurate determination of the average density of baryons in the Universe. 

First, it was shown that there is no plausible astrophysical site for the production of 



deuterium |y], |82|1 ; due to its fragility, post big-bang processes only destroy it. Thus, the 
presently observed deuterium abundance serves as a lower limit to the big-bang production. 
This argument, together with the strong dependence of big-bang deuterium production on 
the baryon density, led to the realization that D is an excellent "baryometer" |T!], Q, and 
early measurements of the deuterium abundance |]84] , B5] , a few parts in 10 5 relative to 



hydrogen, established that baryons could not contribute more than about 20% of closure 
density. This important conclusion still holds today. 



The chemical evolution of 3 He is more complicated. Helium-3 is produced in stars as they 
burn their primeval D before reaching the main sequence, and later by the nuclear reactions 
that cook hydrogen into helium. Some massive stars destroy (or astrate) 3 He. It wasn't until 
the late 1970's that a suitable argument for 3 He was found: The present sum of D + 3 He 
bounds their combined big-bang production [p6| . 

Lithium was the last to come into the fold. Stellar processes both destroy and produce 7 Li; 
moreover, the abundance of 7 Li varies greatly, from 7 Li/H ~ 10 -9 in the interstellar medium 
(ISM) to less than 7 Li/H ~ 10~ 12 in some stars. In 1982, Spite and Spite circumvented these 
difficulties by measuring the 7 Li abundance in the oldest stars in our galaxy, metal-poor, pop 
II halo stars. They found 7 Li/H~ 10~ 10 [JS7H , which is consistent with big-bang production. 



Their results established the case for the primeval 7 Li abundance, which since has been 
strengthened by the work of others |88|, |89| . 



For the last decade much effort has been devoted to the critical comparison of the theo- 
retical predictions and inferred primordial abundances of the light elements. The predictions 
depend upon the ratio of baryons to photons (= rf). As we shall discuss in more detail, if 
r\ is between about 2.5 x 1CT 10 and 6 x 1CT 10 there is concordance between the predicted 
and measured abundances of all four light elements (see Fig. 1). This leads to the best 
determination of the baryon density, 

p B = rin^mN = 1.7 x 1(T 31 gcnT 3 -4.1 x l(r 31 gcm" 3 , (1) 

where the number density of photons, n 7 = 411cm -3 , is known very precisely because the 
CBR temperature is so well determined [p0| , T = 2.726 K ± 0.005 K. On the other hand, 
because the critical density, 

p crit = 3H$/8irG = 1.88/i 2 x 10~ 29 gcnr 3 , (2) 

depends upon the Hubble constant, which is still only known to within a factor of two, the 
fraction of critical density contributed by baryons is less well known: 

tt B = 0.009/r 2 -0.02/i" 2 , (3) 

where h = ifo/100kms _1 Mpc -1 . For a generous range for the Hubble constant, h = 0.4 — 1, 
baryons contribute between 1% and 14% of closure density. 

This fact has two profound implications. First, since "optically" luminous matter (stars 
and associated material) contributes much less than 1% of the critical density, fi LUM ~ 
0.003/i _1 |JT], most baryons must be dark, e.g., in the form of hot, diffuse gas, or "dark 



stars" that have either exhausted their nuclear fuels (black holes, neutron stars or white 
dwarfs) or were not massive enough (less than about O.O8M ) to ignite them. In clusters 
of galaxies most of the baryonic matter seems to be in the form of hot, x-ray emitting gas. 
Further, there is now indirect evidence for the existence of dark stars, known as MACHOs 
for Massive Astrophysical Compact Halo Objects, through their gravitational microlensing 
of distant stars |J2] . 



Second, there is strong — though not yet conclusive — evidence that the average mass 
density of the Universe is significantly greater than 14% of the critical density [j93| ; if this is 



indeed the case, most of the mass density of the Universe must be "nonbaryonic," with the 



most promising possibility being elementary particles left over from the earliest moments of 
the Universe ||94|| . Large-scale experiments are underway in laboratories all over the world 



to directly detect the nonbaryonic dark matter associated with the halo of our own galaxy 



Big-bang nucleosynthesis plays the central role in defining both dark-matter problems 
which are central to cosmology today. For example, the detection of temperature variations 
in the CBR by the COBE satellite was a dramatic confirmation of the general picture that 
structure evolved from small density inhomogeneities amplified by gravity. One of the great 
challenges in cosmology is to formulate a coherent and detailed picture of the formation of 
structure (i.e., galaxies, clusters of galaxies, superclusters, voids, and so on) in the Universe; 
the nature of the dark matter is crucial to doing so. 

Primordial nucleosynthesis also allows us to "study" conditions in the early Universe, and 
thereby, to probe fundamental physics in regimes that are beyond the reach of terrestrial 
laboratories. For example, more than ten years ago the overproduction of 4 He was used to 
rule out the existence of more than three light (mass less than about 1 MeV) neutrino species 
and constrain the existence of other light particle species [^Bl ^, [57] , 



The remainder of this article is given to a careful assessment of the predictions and 
observations, paying special attention to the conclusions that can be sensibly draw about 
the baryon density. We begin with the easier part, a discussion of the theoretical predictions, 
where the few uncertainties are primarily statistical in nature and easy to quantify. We then 
move on to the more difficult side, observations. Here the situation is just the reverse: 
The uncertainties are dominated by possible systematic errors and interpretational issues, 
which cannot be characterized by standard Gaussian error flags; care and judgment must be 
exercised to reach reliable conclusions. 

2 Theoretical expectations 

Things have come a long way from the early semi-quantitative estimates of Gamow and 
his collaborators []?JJ. The first serious attempts to calculate big-bang 4 He production were 
those of Alpher, Follin and Herman |99j and Hoyle and Tayler [76|. After the discovery of the 



CBR, Peebles JT8) , and then Wagoner, Fowler, and Hoyle |7j| , carried out detailed numerical 
calculations incorporating a reaction network of more than ten nuclides (see Fig. 2). The code 
developed by Wagoner et al. [fn|, |100|1 (using numerical techniques introduced by Truran ||101 | 



to study nucleosynthesis in exploding stars) has since been modified to incorporate small, 
but significant corrections (finite-temperature effects, finite nucleon-mass effects, the slight 
coupling of neutrinos to the electromagnetic plasma, and so on), to update nuclear- reaction 
rates, and to improve the numerical-integration scheme. Its present-day incarnation has 
become the standard code (though other codes have been developed independently), and is 
described in detail in Refs. ||102| , |103| . 



The assumptions underlying the standard scenario of big-bang nucleosynthesis are few: 
(i) Big-bang (i.e., Friedmann- Robertson- Walker) cosmological model; (2) three massless (or 
very light) neutrino species; (3) small or vanishing neutrino chemical potentials; (4) no 
additional light particle species present in thermal abundance; (5) spatially uniform baryon- 
to-photon ratio. In addition, there are "nuclear input parameters:" neutron mean lifetime, 



which determines the matrix element for all the reactions that interconvert neutrons and 
protons, and cross sections for nuclear reactions. 

As recently as ten years ago the uncertainty in the mean neutron lifetime was significant; 



at present it is known very precisely: r n = 889 sec ± 2 sec [ |104|| . The other cross sections that 
are required have been measured in the laboratory at energies appropriate for primordial nu- 
cleosynthesis (this is in contrast to stellar nucleosynthesis where lab-measured cross sections 
must be extrapolated to much lower energies). With the exception of 7 Li, the uncertainties 
in cross sections do not result in significant uncertainties in the light-element yields. 

Two recent Monte-Carlo studies have quantified the uncertainties in the predicted abun- 
dances ||103| , |105| . Kernan and Krauss |105|| ran a suite of 1000 models with input parameters 



chosen from the probability distributions for the various cross sections and neutron mean 
lifetime. For a baryon-to-photon ratio of 3 x 10~ 10 the "two-sigma" range of the abundances 
found was: Y P = 0.239-0.241; D/H = 6.7 x 10~ 5 -9.0 x 10" 5 ; 3 He/H = 1.4 x 10" 5 -1.9 x 10~ 5 ; 
and 7 Li/H = 0.81 x 10~ 10 -1.7 x 10~ 10 (i.e., 950 of the models had abundances within the 
stated intervals). Here, Y P is the mass fraction of 4 He produced; the other abundances are 
specified by number relative to hydrogen. Only for 7 Li is the uncertainty significant when 
compared to the uncertainty in the observed abundance. It arises due to three cross sections 
that are still poorly known: 3 He+ 4 He — > 7+ 7 Be, 3 H+ 4 He — > 7+ 7 Li, and 7 Li+p — ► 4 He+ 4 He. 
(In fact, there could be systematic errors in one or more of these cross sections, as the results 
of different experiments are not consistent with their quoted error flags [|103| , |105|| .) 



As an aside, the 4 He yield is known most accurately, apparently to a precision of better 

than 1%. For this reason even tiny corrections have become important, and it is difficult to 

judge whether or not every significant effect has been taken into account. The most recent 

correction may serve as a guide: finite nucleon-mass effects led to an increase, AYp ~ +0.001 

10(j|| . Further, an informal poll of the various nucleosynthesis codes known to us gave results 



that, with the exception of 4 He, fell within the above mentioned Monte-Carlo range; for 4 He, 
the values reported ranged from Yp = 0.237 to Yp = 0.241. 

Modifications of the standard scenario have also been investigated, including almost 
every imaginable possibility |107j| : additional light particle species; unstable, massive tau 



neutrino; decaying particles; variations in the fundamental constants; large neutrino chemical 
potentials; primeval magnetic fields; and spatial variations in the baryon-to-photon ratio. In 
most instances the "nonstandard physics" was introduced for the purpose of obtaining a 
bound or limit based upon primordial nucleosynthesis, e.g., the previous mentioned limit to 
the number of light neutrino species. In a few cases, however, it was motivated by other 
considerations. 

For example, Witten suggested that if the transition from quark/gluon plasma, which ex- 
isted prior to about 10 -5 sec, to hadron matter involved a strongly first-order phase transition 
the resulting distribution of baryons could be quite inhomogeneous ||108|| , thereby possibly sig- 
nificantly changing the outcome of primordial nucleosynthesis. For a short time, it appeared 
that such inhomogeneity could lead to a relaxation of the bound to Qp, even permitting 
closure density in baryons |109| . It is now clear that any significant level of inhomogeneity 



upsets the agreement of the predictions with the observations ||110j , |111| , |112| ; moreover 



there is now little motivation from particle physics for a strongly first-order quark/hadron 
phase transition. At present, the only modification involving the known particles that leads 
to significant changes is the possibility that the tau neutrino has a mass of the order of 



!MeV-30MeV |113|| . The present laboratory limit to its mass is just above 30MeV and 



should be improved enough to clarify this issue soon. 

To summarize the theoretical situation; the only compelling scenario for primordial nu- 
cleosynthesis is the standard one. Within the standard picture the predictions for the light- 
element abundances have uncertainties that, with the exception of 7 Li, are not significant 
when compared to the accuracy with which the primeval abundances are known. At present, 
the comparison between theory and observation turns primarily on the observations, and 
there, the uncertainties are more difficult to quantify. 

3 Confrontation Between Theory and Observation 

The predictions of the standard scenario, including "two-sigma" uncertainties based upon 
our Monte-Carlo calculations are shown in Fig. lfj We now discuss the inferred primordial 
abundances, with emphasis on inferred, since one must deduce the primordial abundances 
of D, 3 He, 4 He, 7 Li from material that has undergone some 10 Gyr of chemical evolution. 

3.1 Deuterium and Helium-3 

Since deuterium is the most weakly bound, stable nucleus it is easy to destroy and difficult 
to produce. As discussed earlier, the deuterium abundance observed today provides a lower 
limit to the big-bang production. The Apollo Solar Wind Composition experiment, which 
captured solar-wind particles in foils exposed on the moon, and the subsequent analysis by 
Geiss and Reeves |84j] provided the first accurate assessment of the pre-solar D and 3 He 



abundances. Based on these experiments and studies of primitive meteorites, Geiss deduces 



a pre-solar (i.e., at the time of the formation of the solar system) deuterium abundance ||114 

D 



=2.6 ±1.0 x l(T b . (4) 

H/ 

This value is consistent with measurements of the deuterium abundance in the local ISM 



(i.e., within a few hundred pc) made two decades ago by the Copernicus satellite [pq] , and 



more recently by the Hubble Space Telescope JTT 

'§) = 1.65185 x 10- 5 . (5) 

n/HST 

That the ISM value is slightly lower than the pre-solar abundance is consistent with slow 
depletion of deuterium with time since the material in the ISM is about 5 Gyr younger than 
the material from which our solar system was assembled. A sensible lower bound to the 
primordial deuterium abundance, 

^) p >l-6xl0" 5 , (6) 



^ur Monte-Carlo calculations are similar to those in Refs. [103, 105] with one exception; we have treated 
the cross sections for 7 Li production differently. For the cross sections where data from two experiments are 
inconsistent we have used both distributions, alternating between the two. 



based on these measurements leads to an upper limit to 77 of 9 x 1CT 10 . Because of the rapid 
variation of the amount of deuterium produced with 77, this upper limit is rather insensitive 
to the exact lower bound adopted for D/H. Further, this argument is very robust because 
it involves minimal assumptions about galactic chemical evolution, simply that D is only 
destroyed by stellar processing |82[ . 



It would be nice to exploit the rapid variation of deuterium production to obtain a lower 
bound to 77, based upon the overproduction of deuterium. This cannot be done directly 
because deuterium is so easily destroyed. However, an equally useful bound can be derived 
based upon the sum of D + 3 He production. Primordial deuterium either resides in the ISM 
or has been burnt to 3 He (via D + p — > 7 + 3 He). A significant fraction of 3 He survives 
stellar processing (in fact, low mass stars are net producers of 3 He), and thus an upper bound 
to the primordial D + 3 He abundance can be inferred from present-day measurements with 



few assumptions. This line of reasoning was introduced by Yang et al. |86| who derived the 
bound, 



H Jp V H /© V H /© 



where the He survival fraction g 3 was argued to be greater than 25% |116|. (It should be 
noted that even massive stars, which tend to burn 3 He, eject some 3 He in their winds.) The 
bound was improved by taking account of material that has been processed by more than 



one generation of stars |97|, |98| . Both methods lead to similar upper limits to the primordial 
D + 3 He abundance, 

/D + 3 He\ . 

HH,- 11 * 10 ■ (8) 

which in turn gives the bound 77 > 2.5 x 10~ 10 . Like the upper limit to rj based upon 
deuterium, this lower limit is insensitive to the precise bound to the primeval abundance of 
D + 3 He because of the steep rise of D + 3 He production with decreasing 77. Together, D 
and 3 He define a concordance interval, 2.5 x 10~ 10 < 77 < 9 x 10 -10 . 

The theoretical belief that low-mass stars actually increase the D + 3 He abundance by 
producing 3 He is supported by the observations of Wilson, Rood and Bania [ p.!7| |. By using 
the analogue of the 21 cm hydrogen hyperfine transition for 3 He they found 3 He/H ~ 10~ 3 
in planetary nebulae. This much additional 3 He production agrees with the value predicted 
by stellar models of Iben and Truran [118] - However, measurements of the 3 He abundance 
by the same method in hot, ionized gas clouds, so called H II regions, vary greatly, from 
3 He/H ~ 1 x 1CT 5 to 3 He/H ~ 8 x 1CT 5 ||119| , which suggests that 3 He is destroyed by 



varying degrees |120| . While H II regions are the only place outside the solar system where 
the 3 He abundance can be measured, they are samples of the cosmos dominated by the 
effects of massive, young stars, which are the most efficient destroyers of 3 He, and thus, they 
do not represent "typical samples" of the cosmos (so far as the chemical evolution of 3 He is 
concerned). In any case, we believe that a 3 He survival fraction of 25% or more remains a 
reasonably conservative estimate as applied to the solar system 3 He abundance. 



3.2 Helium-4 

In two important regards the primordial 4 He abundance is the easiest to measure: it is large, 
around 25% by mass fraction, and the chemical evolution of 4 He is straightforward — stars are 
net producers of 4 He. On the other hand, the predicted abundance is most accurately known 
and varies only logarithmically with rj. Thus, measuring the 4 He abundance to sufficient 
accuracy to sharply test the big-bang prediction is just as challenging as determining the 
other light-element abundances. 

Needless to say, since 4 He is ubiquitous, its abundance can be measured in many different 
ways, all of which give values consistent with a primeval mass fraction of around 25%. The 
most accurate determinations of the primeval 4 He abundance rely on measurements of its re- 
combination radiation in low-metallicity H II regions (see Ref. [|121|| for a detailed discussion 



of the experimental method). Since stars produce both helium and other heavier elements, 
contamination due to stellar production should be minimized in metal-poor samples of the 
Universe. A number of groups have obtained high-quality data for very metal-poor, extra- 
galactic H II regions, which has allowed determination of helium abundances to very good 



statistical accuracy [122]. Moreover, several independent and detailed analyses of these data 



sets have been carried out PSL |123|, |124|. The quality of the data and the accuracy of the 



abundance determinations desired are now such that possible systematic errors dominate the 
error budget, and they are the focus of our attention. 

3.2.1 Systematic Effects 

The first step in the path to the primordial helium abundance is measuring line strengths of 
the recombination radiation for hydrogen and helium. Line strengths are then translated into 
a helium mass fraction by means of theoretical emissivities for both helium and hydrogen 
and modeling of the H II region. In modeling an H II region spherical symmetry and uniform 
temperature are assumed, neither of which is an excellent assumption given that a typical 
H II region is heated by a few massive, young stars near its center. Since the ionization 
potentials for hydrogen and helium are different, corrections must be made for both neutral 
and doubly ionized helium. Collisional excitation can be significant, but is not easy to 
accurately estimate. Stellar absorption by the stars heating the H II region can affect the 
excitation of the hydrogen and helium in the H II region. Absorption due to intervening dust 
can also affect the abundance determinations. A summary of our estimate of the systematics, 
based largely on the discussion of Skillman and Kennicutt [|125|| and Skillman et al. ||126 



is given in Table |l|. A detailed numerical assessment of some of these effects has recently 
been carried out by Sasselov and Goldwirth ||127|| , who suggest that the systematic errors 
may even be slightly larger. 

3.2.2 Primordial helium-4 abundance 

Even in the most metal-poor H II regions some of the 4 He is produced by stars. Since stars 
also produce the elements beyond 4 He (collectively referred to as metals), there should be 
a direct relationship between metallicity and stellar-produced 4 He. Peimbert and Torres- 
Peimbert [|128|| pioneered the extrapolation of the helium abundance vs. heavy-element 



abundance to zero metallicity to infer the primordial 4 He abundance. Oxygen, nitrogen, and 
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Type of correction Estimate 

Line ratios (including dust absorption) ±2% 

Emissivities ±2% 

Collisional excitation and stellar absorption ±1% 

Neutral helium +2% 

"Total +7%,-5% 

Table 1: Estimate of systematic errors. 

carbon have all been used as indicators of stellar nucleosynthesis and hence the amount of 
stellar produced 4 He. Each has its advantages and disadvantages ||129|] , though the quanti- 
tative results are very similar. Recently, Olive and Steigman |124j| have performed a detailed 



statistical analysis of very metal-poor H II regions, and derive a primordial 4 He abundance 
(see Fig. 3) 

Y P = 0.232 ± 0.003lg;gi2 ( 9 ) 

where their statistical error is quoted first and the systematic error based upon Table [I] 
appears second. (For reference, their quoted systematic error is ±0.005.) 

To summarize, there is undisputed evidence for a large primeval 4 He abundance whose 



only plausible explanation is the big bang. Following Olive and Steigman ||124|| we take as 
a reasonable estimate for the primeval mass fraction, Yp = 0.221 — 0.243, which allows for 
a two-sigma statistical uncertainty and their one-sigma systematic uncertainty. Within the 
two-sigma theoretical uncertainty, such a primeval mass fraction of 4 He is consistent with the 
big-bang prediction provided 0.8 x 10~ 10 < r] < 4 x 10~ 10 . At present, errors are dominated 
by possible systematic effects; allowing for our higher estimate of systematic error, a primeval 
4 He mass fraction as low as 0.214, or as high as 0.254, cannot be ruled out with certainty. 
This extreme range for the primeval 4 He abundance is consistent with a much larger interval, 
6 x 10~ n < r] < 1.5 x 10 -9 , illustrating the exponential sensitivity of 7] to Yp. 

3.3 Lithium 

The study of extremely metal-poor, pop II halo stars has provided the bulk of our knowl- 
edge of the light elements beyond 4 He. It began with the work of Spite and Spite |j57| , 
who measured the 7 Li abundance as a function of metallicity (iron abundance) and sur- 
face temperature. Much to their surprise they found a "plateau" in the 7 Li abundance and 
established what has become a very strong case for the determination of the primeval 7 Li 
abundance. 

The Spite plateau refers to the fact that the 7 Li abundance as a function of surface 
temperature is flat for surface temperatures greater than about 5600K (see Fig. 4), and 
further that the 7 Li abundance for stars with temperatures greater than 5600K as a function 
of iron abundance is flat for very low iron abundance (Fig. 5). The first plateau gives 
strong evidence that the stars with the highest surface temperatures are not destroying 
their 7 Li by convective burning since the depth of the convective zone depends upon the 
surface temperature (for temperatures lower than 5600K the measured 7 Li does vary with 



surface temperature indicating convective burning). The second plateau indicates that any 
7 Li due to stellar production must be insignificant for the most metal-poor stars since the 
7 Li abundance is independent of the metal abundance. 

The actual value of the 7 Li abundance on "the Spite plateau" is subject to several impor- 
tant systematic effects. In particular, model atmospheres used by different authors assume 
effective surface temperatures, differing by as much as 200 K. Other differences in the model 
atmospheres, including assumptions made about opacities, also affect the inferred 7 Li abun- 
dances in a systematic way. These systematic effects explain the main difference between 
the Spite and Spite abundance, 7 Li/H = 1.1 x 10 -10 , and the value derived recently by 
Thorburn p9| from a sample of 90 pop II stars, 7 Li/H = 1.7 x 10~ 10 (see Figs. 4, 5). Fur- 



ther, Thorburn's data seems to indicate a slight systematic variation of the 7 Li abundance 
with surface temperature, possibly indicating some depletion from a higher primordial value 
by processes that transport 7 Li inward to regions of high enough temperature that it can 
be burned; e.g., meridonal mixing |130|. However, the amount of depletion is constrained 



by the relatively narrow spread in 7 Li abundance for a wide range of surface temperatures 
and metallicities. Microscopic diffusion is ruled out by this fact, though stellar models that 
incorporate rotation, which suppresses microscopic diffusion, can be made consistent with 



the observations 131 



The case against significant depletion (and hence for a primeval abundance) was further 
strengthened by the observation of 6 Li in a pop II star by Smith, Lambert, and Nissen [|132[ ; 



Hobbs and Thorburn [ |133j| have detected Li in this and another pop II star. Big-bang 
production of 6 Li is negligible and so the 6 Li seen was presumably produced by cosmic-ray 
processes, along with beryllium and boron (as discussed below). Since 6 Li is much more 
fragile than 7 Li and yet still survived with an abundance relative to Be and B expected for 
cosmic-ray production, depletion of 7 Li cannot have been very significant ||134|, |135||. These 



6 Li measurements allow for a largely model-independent limit to the amount of 7 Li depletion, 
less than about a factor of two. 

In summary, based on metal-poor, pop II halo stars we infer a primordial 7 Li abundance 
of 

(j?) =1.4±0.3^jxHT 10 (10) 

where the central value is the average of the Spite and Spite and Thorburn determinations, 
the statistical error is listed first, and the systematic error second. The systematic-error 
estimate consists of ±0.4 due to differences in model atmospheres and +1.4 to account for 
possible depletion. In fixing a range for the primordial 7 Li abundance it is the systematic 
error that is most important; accordingly, we use the sum of statistical plus systematic error 
to derive our estimate for the 7 Li abundance, 0.7 x 10~ 10 < 7 Li/H < 3.5 x 10 -10 . Allowing for 
the two-sigma theoretical uncertainty, the concordance interval is 1 x 10~ 10 < r] < 6 x 10~ 10 . 
(We note that the 95% confidence range for the 7 Li abundance advocated by Thorburn [|5T5| 
differs only slightly from ours.) 

3.4 Beryllium and Boron 

While the inhomogeneous variant of big-bang nucleosynthesis motivated by a first-order 
quark/hadron phase transition cannot alter the basic conclusions, an important question 

9 



remains, namely, is there an observable signature that can differentiate between the in- 
homogeneous and the homogeneous models, thereby probing the quark/hadron transition? 



Several authors ||109| argued that the regions with high neutron-to-proton ratio that exist in 
inhomogeneous models could lead to "leakage" beyond mass 5 and mass 8 and traces of 9 Be, 
10 B, n B, and possibly even r-process elements (neutron-rich isotopes) could be produced. 
However, detailed studies by Sato and Terasawa [ |1 1 1|| and Thomas et al. ||112| have shown 



that such leakage is negligible when the D, 3 He, 4 He, and 7 Li abundances are consistent with 
their observed values, with Be/H, B/H ~ 10 -18 

An observational program similar to that of Spite and Spite |87|| was begun for beryllium 



and boron. Recently, both beryllium [|136|| and boron ||137|| have been detected in metal-poor 



pop II halo stars. The observations indicate that beryllium and boron abundances scale 
with metallicity, strongly suggesting that their production was not the big bang ||138| . The 
processes that produce the beryllium and boron (and 6 Li) seen in younger pop I stars (like 
our sun) are thought to be cosmic-ray reactions | |139| |. For Be and B, such reactions involve 
the breakup of heavy nuclei such as C, N, O, Ne, Mg, Si, S, Ca, and Fe by protons and alpha 
particles (for lithium in pop II stars, alpha plus alpha fusion reactions are dominant |140| |). 



3.5 Toward truly primordial abundances 

As the reader by now should appreciate, the task of disentangling 10 Gyr of galactic chemical 
evolution is not an easy one. What are the prospects for determining the light-element 
abundances in very primitive samples of the Universe (that is, in objects seen at very high 
red shift)? 

Hydrogen clouds at high redshift "backlit" by quasars offer the possibility of measur- 
ing the deuterium abundance in very old, very distance, and very primitive samples of the 
|1 11|| . These clouds, known as quasar absorption line systems, are "seen" by the 



cosmos 



absorption features they produce in the quasar spectrum; many are observed to be very 
metal-poor, indicative of primeval material. There have been many searchs for the deu- 
terium analog of the 1216 A Lyman-a absorption feature, which is shifted very slightly to 
the blue, by about 0.33 A. Recently, Songaila et al. and Carswell et al. ||142| announced a 



possible detection of deuterium in a redshift z = 3.32 absorption line system (in the quasar 
Q00 14+813); if it is deuterium, it corresponds to an abundance 

~) = (1.9-2.5) x lO" 4 . (11) 

Both groups are quick to point out that a single measurement does not constitute a definite 
detection of deuterium as there is a significant probability (~ 15%) that the feature seen 
arises from Lyman-a absorption due to another, smaller hydrogen cloud at slightly lower 
redshift. At the very least though, their detection provides an important upper limit to the 
deuterium abundance in this very primitive sample of the cosmos. 

Interpreting the detection as an upper bound to the primordial deuterium abundance 
leads to the constraint rj > 1.6 x 10~ 10 , only slightly less stringent than the previous bound 
based upon the production of D + 3 He. If, on the other hand, it is interpreted as a measure- 
ment of the primordial deuterium abundance, then (D/H) p ~2x 10~ 4 ^> [(D + 3 He)/HL ~ 
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4 x 10~ 5 , which leads to a problem in understanding the observed pre-solar D+ 3 He abun- 
dance. Since it is almost certain that D is destroyed by burning to 3 He one would expect a 
much higher D + 3 He abundance than has been observed. This could indicate a problem with 



models of the chemical evolution of He, or the interpretation as a deuterium detection ||143 

Carswell et al. have studied another quasar at a similar redshift (Q0420-388) and have 
detected deuterium in an absorption line system at the level of about 2 x 10 -5 , with a three- 
sigma upper limit of 6 x 10~ 5 |144f| . This new observation seems to imply that the previous 
detection was really due to another small hydrogen cloud. Further, if correct, it fits nicely 
with the local measurements of D and 3 He, and suggests that the astration of D is not great. 

The merits of using quasar absorption line systems to obtain the primordial deuterium 
abundance are clear. The means necessary for such observations are now at hand: large- 
aperture telescopes with very high-resolution spectrometers, such as the 10 meter Keck tele- 
scope used by Songalia et al. With some luck, it should just be a matter of time before 
the deuterium Lyman-a; feature is measured in several such systems. Once it is, and if the 
abundances are similar, both the reality of the feature and its interpretation as reflecting the 
primordial D abundance will have been established. 

With regard to 3 He, one might hope to eventually use the 3 He hyperfine line to determine 
the 3 He abundance in extragalactic H II regions that are very metal poor. However present 
technology is only marginally sufficient to observe galactic 3 He so it will take time before 
extragalactic detections are possible. 

The 4 He abundance has been measured through its absorption lines in a quasar at redshift 
: = 2.72 (HS1700+6414) ||145|| , and, very recently, observations made with the refurbished 



Hubble Space Telescope have revealed the presence of singly ionized 4 He in the intergalactic 
medium ||146| . While both measurements provide important confirmation of a large, primeval 



4 He abundance in very primitive samples of the cosmos, they lack the precision necessary to 
sharply test big-bang nucleosynthesis. In that regard, metal-poor extragalactic H II regions 
provide the most accurate determinations. 

Finally, owing to its small abundance, it seems very unlikely that the 7 Li abundance can 
be measured in high-redshift objects, or even in extragalactic stars. On the other hand, the 
data at hand present a good case for having determined the 7 Li abundance in the very oldest 
stars in our galaxy. 

4 Implications and Future Directions 

The agreement between the predictions of primordial nucleosynthesis and the inferred pri- 
mordial abundances is impressive — all the more so when viewed in light of the sharpening 
of the theoretical predictions and the improvement in the observational data that the past 
decade has witnessed. Without a doubt, primordial nucleosynthesis provides the most sig- 
nificant test of the standard cosmology — which it passes with flying colors — and leads to the 
best determination of the density of ordinary matter. 

Where do we stand? The data are not yet good enough to single out a value for the 
baryon-to-photon ratio. However, they are good enough to delineate a very narrow "con- 
cordance interval" where the predicted abundances of all four light elements are consistent 
with their measured values. 



11 



The lower limit to the concordance interval hinges primarily upon the D + 3 He abundance. 
Based upon our understanding of the difficulty of efficiently destroying 3 He, 77 = 2.5 x 1CT 10 
stands as a reliable lower bound. This lower bound is buttressed by both 7 Li — for 77 < 
1 x 10~ 10 the predicted 7 Li abundance rises above 3.5 x 1CT 10 — and by the Songaila et 
al. upper limit to the primitive deuterium abundance — for 77 < 1.6 x 1CT 10 D/H exceeds 

2.5 x io- 4 .g 

The upper limit to the concordance interval derives from 4 He, 7 Li and D, with the strin- 
gency of the limits in that order, but the reliability in the reverse order. Assuming that the 
primordial mass fraction of 4 He is no larger than 0.243, based upon the work of Olive and 
Steigman [|124j| , then 77 must be less than 4 x 10~ 10 . On the other hand, if owing to systematic 



error Yp is as large as 0.254, then 77 could be as large as 1.5 x 10 9 . This illustrates the point 
mentioned earlier, the upper limit to 77 depends exponentially upon the upper limit to Yp, 

77 max -4x 10~ 10 exp[100(Y^ ax - 0.243)]. (12) 

While 4 He is arguably the most striking confirmation of big-bang nucleosynthesis, the loga- 
rithmic dependence of the 4 He mass fraction on 77 makes it a very poor baryometer. 

The uncertainty in our reasonable upper bound to 7 Li, 7 Li/H < 3.5 x 10~ 10 , is primarily 
systematic error associated with possible 7 Li depletion in metal-poor, pop II stars. Our 
upper bound to 7 Li implies 77 < 6 x 10~ 10 . On the other hand, since the strongest argument 
against very significant depletion of 7 Li in metal-poor, pop II stars is the observation of 6 Li, 
which has only been seen in two stars, very significant depletion of 7 Li cannot be ruled out 
with certainty. Taking as an extreme upper limit, 7 Li/H~ 6 x 10~ 10 , corresponding to a 
factor of four depletion, 77 could be as large as 9 x 10~ 10 . 

Finally, turning to deuterium; because it is so easily destroyed and lacks a plausible 
contemporary astrophysical site for its origin, its primordial abundance must be larger than 
what is seen today: D/H > 1.6 x 10~ 5 . This implies an upper bound to 77 of 9 x 10~ 10 . It 
seems very difficult to get around this simple argument; moreover, because the abundance 
of deuterium varies so rapidly with 77 this bound is insensitive to the precise deuterium 
abundance assumed. 

The difficulty in drawing sharp conclusions from the comparison between predicted and 
measured primordial abundances — e.g., stating two- or three-sigma limits — is the fact that 
the dominant uncertainties, primarily in the observations, are not Gaussian statistical errors. 
Two- or three-sigma limits in the present circumstance simply have no meaning. Instead, we 
choose to quote a "sensible" and an "extreme" concordance interval for the baryon-to-photon 
ratio. For the sensible interval we take 2.5 x 10~ 10 to 6 x 10~ 10 , supported from below by 
D+ 3 He overproduction and above by 7 Li overproduction. Some have argued that 4 He can 
be used to push the upper limit down to 4 x 10~ 10 ; however, as described, such an upper 
limit is exponentially sensitive to the uncertainties in the primeval 4 He abundance, and thus 
not very robust. 

It is interesting to compare our sensible interval with similar concordance intervals found 



by Yang et al. [86 in 1984, 4 x 10" 1U < 77 < 7 x 10~ iU , and by Walker et al. |8| in 1991 



2.8 x 10 < 77 < 4 x 10 . The difference between our lower limit and that of Yang et 



2 Had we used the Carswell et al. [144] upper limit, D/H < 6 x 10 5 , which has yet to be published, the 
lower bound to r\ would be very similar to that based upon D + 3 He. 
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al. is the fact that we have allowed for slightly more astration of 3 He. The somewhat larger 
difference between our upper limit and that of Walker et al. involves 7 Li: They used the 
lower Spite and Spite value for the primordial 7 Li abundance and did not allow for systematic 
error. In any case, the fact that the differences between the concordance intervals are small 
is reassuring. 

In setting the extreme range, we take account of of our less than perfect understanding 
of the chemical evolution of the Universe during the 10 Gyr or so since primordial nucle- 
osynthesis, as well as other possible systematic errors. Though there is no plausible reason 
for believing that 3 He could be astrated significantly, or that the primeval 7 Li abundance is 
significantly different from that seen in halo pop II stars, we do not believe the wider interval 
of r\ = 1.6 x 10~ 10 to 9 x 10~ 10 can be excluded with absolute certainty. Our extreme range 
derives exclusively from deuterium: the present abundance, D/H> 1.6 x 10~ 5 , and the limit 
to the primeval abundance, D/H< 2.5 x 10~ 4 . Moreover, for such extreme values of r] all 
the light-element abundances are pushed to almost untenable values. 

Based upon these concordance intervals — sensible and extreme — we can obtain bounds 
to the baryonic fraction of critical density, albeit at the expense of additional dependence 
upon the Hubble constant, 

sensible : 2.5 x 1(T 10 < 77 < 6 x 1(T 10 => 0.009 < 0.009/T 2 < Q B < 0.02/T 2 < 0.14 
extreme : 1.6 x 10~ 10 < r) < 9 x 10" 10 => 0.006 < 0.006/T 2 < Q B < 0.03/T 2 < 0.21 

where the outer limits to Qb allow for 0.4 < h < 1 (see Fig. 6). 

The implications of these bounds for cosmology are manifold and very significant. First 
and foremost, the nucleosynthesis limit is pivotal to the case for both baryonic and non- 
baryonic dark matter. The nucleosynthesis determination of the baryonic fraction of critical 
density taken together with the observational data that indicate that luminous matter con- 
tributes much less than 1% of critical density and that the total mass density is greater than 
14% of critical density makes the case for these two most pressing problems in cosmology. 
(The Hubble-constant dependence of the nucleosynthesis limits precludes addressing both 
problems by simply choosing the right value for the Hubble constant since both the upper 
and lower limits to Qb scale in the same way.) 

Second, one can exploit the relatively well known baryon density to estimate the total 
mass density by using measurements of the ratio of total mass-to-baryonic mass in clusters of 
galaxies, as determined by recent x-ray studies made using the ROSAT satellite. Assuming 
that rich clusters like Coma provide a fair sample of the "universal mix" of matter (if some, 
or all, of the nonbaryonic dark matter is neutrinos of mass 5 eV to 30 eV, this might not 
be the case) and that the hot, x-ray emitting gas is in virial equilibrium supported against 



gravity only by its thermal motion, White et al. [p. 4711 infer a total mass-to-baryonic mass 



ratio of (20 ± 5) /i 3//2 , which leads to the following estimate for the total mass density: 

^ M TO T(Coma)^ . ^ , i /9 , s 

n ° = ~wih — v Qb - °- 15 - °- 5 h • 13 

Mb (Coma) 

If h is near the lower extreme of current measurements this determination of Q$ lends some 
support to the theoretically attractive notion of a flat Universe (i.e., Qq = 1). 
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The "cluster-inventory" estimate of Qq is a new and potentially very powerful method for 
estimating the mean density of the Universe. There are still important systematic sources 
of error and a key assumption. The key assumption is that the baryons are either in stars 
(visible matter) or hot, x-ray emitting gas (by a wide margin, the baryons in the hot gas 
outweigh those in stars). If there is a large amount of baryonic matter hidden in dark stars, 
then Mtot/^b would be smaller. On the other hand, essentially all systematic sources of 
error go in the direction of increasing M^ot/Mb- For example, if the hot gas is partially 
supported by magnetic fields or bulk motion of the gas, then Mtot would be larger. There 
is some evidence that Mtot has been underestimated: The measurement of the mass of 
one cluster of galaxies based upon weak-gravitational lensing of galaxies behind the cluster 
yields a mass than is almost a factor of three larger than that based upon x-ray studies 
||148| . If the hot gas is clumpy, rather than smooth as is assumed, then the gas mass would 



be smaller, which also increases M^q^/Mb- It is intriguing that a factor of two or three 
increase in Mt OT /Mb would bring the estimate of f2 close to unity. In any case, further 
study of the x-ray and gravitational lensing data as well as the better x-ray temperatures 
that the ASCA satellite is now providing should help refine this new technique for estimating 
the mean density. 

While the primary concern of this paper is the baryon density of the Universe, big-bang 
nucleosynthesis also places an important constraint to the number of light particle species 
present around the time of nucleosynthesis, usually quantified as a limit to the equivalent 
number of neutrino species, N v . This limit arises because more species lead to additional 
4 He production [Q. The limit to N u relies upon a lower limit to i] and an upper limit to 
Yp. Using the D + 3 He bound to rj, rj > 2.5 x 10~ 10 , and our reasonable upper limit to 
Y P < 0.243, it follows that N v < 3.4.0 This limit depends upon the upper limit to Y P . 
However, in contrast to the upper limit to rj based upon Y P , the dependence is linear, not 
exponential, N v < 3.4 + {Yp' AX — 0.243)/0.012, and so the limit to N v is far less sensitive to 
Yp than is upper bound to r/. 

Finally, what does the future hold? We believe that primordial nucleosynthesis is the 
best method for determining the mean baryon density, and, in that regard, that deuterium 
is the best baryometer. Not only does the primeval deuterium abundance vary rapidly with 
the baryon-to-photon ratio (as r/ -15 in the relevant range), but prospects for measuring 
its primeval abundance in high-redshift, metal-poor quasar absorption line systems look 
promising. While the present situation is unsettled, with a reported detection of D/H ~ 
2 x 10~ 4 , as well as an upper limit, D/H < 6 x 10~ 5 , the situation should improve. A 
handful of such measurements could establish the primeval D abundance to an accuracy 
of 10%, which would determine the baryon density to better than 5% (taking account of 
both the observational and theoretical uncertainties). Because of their weak dependence 
upon the baryon-to-photon ratio, as well as lingering systematic uncertainties, 4 He and 7 Li 
are destined to play a supporting role, albeit a very important one. It is both ironic and 
satisfying that after twenty years deuterium is still the best baryometer. 

More than forty years have passed since Gamow's introduction of the notion of big-bang 



3 Kernan and Krauss |105| point out that by using the correlations between the theoretical uncertainties 
in 4 He and D + 3 He, one can improve this limit very slightly, by about 0.1 neutrino species, the equivalent 
of reducing Yp by about 0.001. 



14 



nucleosynthesis, and thirty years have passed since the cosmic background radiation was 
discovered. After more than two decades of careful comparison of theory with observation, 
primordial nucleosynthesis has become the earliest and most important test of the standard 
cosmology. Further, it leads to the best measurement of the density of ordinary matter in 
the Universe and provides a powerful laboratory for studying both the early Universe and 
fundamental physics. 
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Figure Captions 

Figure 1: The predictions of big-bang nucleosynthesis. The broken curves indicate the 2a 
theoretical uncertainties based upon our Monte-Carlo analysis. The 4 He abundance is given 
as mass fraction; the other abundances are number relative to hydrogen. The boxes indicate 
the range of baryon-to-photon ratio consistent with the primeval light-element abundances; 
the 4 He box is dotted to remind the reader that 4 He has not been used to derive an upper 
limit to r] because of the exponential dependence of such a limit to Y P (see text). Our 
sensible concordance range, 2.5 x 1CT 10 < r] < 6 x 10~ 10 , comes from D + 3 He and 7 Li 
overproduction. 

Figure 2: The nuclear reaction network used for big-bang nucleosynthesis; the most impor- 
tant reactions are numbered. The broken boxes for mass 5 and 8 indicate that all nuclides 
of this mass are very unstable. 

Figure 3: The 4 He mass fraction vs. nitrogen abundance for very metal-poor, extragalac- 
tic H II regions. The solid line is the best fit to the data from the analysis of Olive and 
Steigman ||134| . (In deriving their fit they did not include some of the higher metallicity H II 



regions) . 

Figure 4: The 7 Li abundance as a function of surface temperature for very metal-poor, pop 
II halo stars. The decreasing 7 Li abundance in the stars with lower surface temperatures 
indicates they have burned some of their 7 Li (consistent with the fact that such stars are 
predicted to have deeper convection zones). The solid and broken lines indicate the Thorburn 
and Spite plateaus respectively. 

Figure 5: The 7 Li abundance as a function of iron abundance (relative to that seen in 
the solar system) for stars with surface temperatures greater than 5600 K. The increase in 
7 Li abundance seen for the stars with higher iron abundance is indicative of additional 7 Li 
due to cosmic-ray processes and stellar production. The solid and broken lines indicate the 
Thorburn and Spite plateaus respectively. For comparison, the abundances of beryllium and 
boron in metal-poor, pop II halo stars are also shown (from Refs. |138|| ). Unlike the 7 Li 



abundance, the B and Be abundances increase with increasing metal abundance, indicative 
of post-big-bang production. 

Figure 6: The fraction of critical density contributed by baryons as a function of the Hubble 
constant for the sensible concordance range of baryon-to-photon ratio (solid) and extreme 
concordance range (dotted). 
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